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Abstract — We  formulate  and  study  the  thinnest  path 
problem  for  secure  communication  in  wireless  ad  hoc 
networks.  The  objective  is  to  find  a  path  from  a  source  to  its 
destination  that  results  in  the  minimum  number  of  nodes 
overhearing  the  message  by  a  judicious  choice  of  relaying 
nodes  and  their  corresponding  transmission  powers.  We 
adopt  a  directed  hypergraph  model  of  the  problem  and 
establish  the  NP-completeness  of  the  problem  in  2-D  net¬ 
works.  We  then  develop  two  polynomial-time  approxima¬ 
tion  algorithms  that  offer  and  -  =  approximation 

ratios  for  general  directed  hypergraphs  (which  can  model 
non-isotropic  signal  propagation  in  space)  and  constant  ap¬ 
proximation  ratios  for  ring  hypergraphs  (which  result  from 
isotropic  signal  propagation).  We  also  consider  the  thinnest 
path  problem  in  1-D  networks  and  1-D  networks  embedded 
in  a  2-D  field  of  eavesdroppers  with  arbitrary  unknown 
locations  (the  so-called  1.5-D  networks).  We  propose  a 
linear-complexity  algorithm  based  on  nested  backward 
induction  that  obtains  the  optimal  solution  for  both  1- 
D  and  1.5-D  networks.  This  algorithm  does  not  require 
the  knowledge  of  eavesdropper  locations  and  achieves  the 
best  performance  offered  by  any  algorithm  that  assumes 
complete  location  information  of  the  eavesdroppers. 

Index  Terms — Hypergraph,  Thinnest  Path,  Secure  Com¬ 
munication,  NP-complete,  Approximation  Algorithms,  Ap¬ 
proximation  Ratio 

I.  Introduction 
A.  The  Thinnest  Path  Problem 

We  consider  the  thinnest  path  problem  for  secure 
communication  in  wireless  ad  hoc  networks.  For  a  given 
source  and  a  destination,  the  thinnest  path  problem  asks 
for  a  path  from  the  source  to  the  destination  that  results 
in  the  minimum  number  of  nodes  hearing  the  message. 
Such  a  path  is  achieved  by  carefully  choosing  a  sequence 
of  relaying  nodes  and  their  corresponding  transmission 
powers. 

At  first  glance,  one  may  wonder  whether  the  thinnest 
path  problem  is  simply  a  shortest  path  problem  with  the 
weight  of  each  hop  given  by  the  number  of  nodes  that 
hear  the  message  in  that  hop.  Realizing  that  a  node  may 
be  within  transmission  range  of  multiple  relaying  nodes 
and  should  not  be  counted  multiple  times  in  the  total 
weight  (referred  to  as  the  width)  of  the  resulting  path, 
we  see  that  the  thinnest  path  problem  does  not  have  a 
simple  cost  function  that  is  summable  over  edges.  But 
rather,  the  width  of  a  path  is  given  by  the  cardinality  of 
the  union  of  all  receiving  nodes  in  each  hop,  which  is  a 
highly  nonlinear  function  of  the  weight  of  each  hop.  One 
may  then  wonder  whether  we  can  redefine  the  weight  of 
each  hop  as  the  number  of  nodes  that  hear  the  message 
for  the  first  time.  Such  a  definition  of  edge  weight  indeed 
leads  to  a  summable  cost  function.  Unfortunately,  in 
this  case,  the  edge  weight  cannot  be  predetermined  until 
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the  thinnest  path  from  the  source  to  the  destination  in 
question  has  already  been  established. 

A  more  fundamental  difference  between  the  thinnest 
path  and  the  shortest  path  problems  is  that  the  thinnest 
path  from  a  single  source  to  all  other  nodes  in  the 
network  do  not  form  a  tree.  In  other  words,  the  thinnest 
path  to  a  node  does  not  necessarily  go  through  the 
thinnest  path  to  any  of  its  neighbors.  The  loss  of  the  tree 
structure  is  one  of  the  main  reasons  that  the  thinnest  path 
problem  is  much  more  complex  than  the  shortest  path 
problem.  Indeed,  as  shown  in  this  paper,  the  thinnest 
path  problem  is  NP-complete,  which  is  in  sharp  contrast 
with  the  polynomial  nature  of  the  shortest  path  problem. 

Another  aspect  that  complicates  the  problem  is  the 
choice  of  the  transmission  power  at  each  node  (within 
a  maximum  value  that  may  vary  across  nodes).  In  this 
case,  the  network  cannot  be  modeled  as  a  simple  graph  in 
which  the  neighbors  of  each  node  are  prefixed.  In  this 
paper,  we  adopt  the  directed  hypergraph  model  which 
easily  captures  the  choice  of  different  neighbor  sets  at 
each  node.  While  a  graph  is  given  by  a  vertex  set  V  and 
an  edge  set  E  consisting  of  cardinality-2  subsets  of  V,  a 
hypergraph  [1]  is  free  of  the  constraint  on  the  cardinality 
of  an  edge.  Specifically,  any  non-empty  subset  of  V  can 
be  an  element  (referred  to  as  a  hyperedge)  of  the  edge  set 
E.  Hypergraphs  can  thus  capture  group  behaviors  and 
higher-dimensional  relationships  in  complex  networks 
that  are  more  than  a  simple  union  of  pairwise  relation¬ 
ships  [2].  In  a  directed  hypergraph  [3],  each  hyperedge  is 
directed,  going  from  a  source  vertex  to  a  non-empty  set 
of  destination  vertices.  An  example  is  given  in  Fig.  l-(a) 
where  we  have  2  directed  hyperedges  rooted  at  a  source 
node  v  with  each  hyperedge  modeling  a  neighbor  set 
of  v  under  a  specific  power.  The  directed  hypergraph 
model  of  the  thinnest  path  problem  is  thus  readily  seen: 
rooted  at  each  node  are  multiple  directed  hyperedges, 
each  corresponding  to  a  distinct  neighbor  set  feasible 
under  the  maximum  transmission  power  of  this  node. 
The  problem  is  then  to  find  a  minimum-width  hyperpath 
from  the  source  to  the  destination  where  the  width  of  a 
hyperpath  is  given  by  the  cardinality  of  the  union  of  the 
hyperedges  on  this  hyperpath. 

B.  Main  Results 

Based  on  the  directed  hypergraph  formulation,  we 
show  that  the  thinnest  path  problem  in  2-D  networks 
is  NP-complete  even  under  a  simple  disk  propagation 
model.  This  result  is  established  through  a  reduction 
from  the  minimum  dominating  set  problem  in  graphs,  a 
classic  NP-complete  problem.  The  most  challenging  part 
of  this  reduction  is  to  show  the  reduced  problem  is  real¬ 
izable  under  a  2-D  disk  model.  We  further  establish  that 
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even  with  a  fixed  transmission  power  at  each  node  (in 
this  case,  the  resulting  hypergraph  degenerates  to  a  stan¬ 
dard  graph),  the  thinnest  path  problem  is  NP-complete. 
We  then  propose  two  polynomial-time  approximation 
algorithms  that  offer  and  -^4 ===  approximation 
ratios  for  general  directed  hypergraphs  (which  can  model 
non-isotropic  signal  propagation  in  space)  and  constant 
approximation  ratios  for  ring  hypergraphs  (which  result 
from  isotropic  signal  propagation).  Here  n  is  the  total 
number  of  vertices. 

We  also  establish  the  polynomial  nature  of  the  prob¬ 
lem  in  1-D  and  1.5-D  networks,where  a  1.5-D  network 
is  a  1-D  network  embedded  in  a  2-D  field  of  eaves¬ 
droppers  with  arbitrary  unknown  locations.  We  propose 
an  algorithm  based  on  a  nested  backward  induction 
(NBI)  starting  at  the  destination.  We  show  that  this  NBI 
algorithm  has  O(n)  time  complexity.  Since  the  size  of 
the  input  data  is  O(n),  the  proposed  algorithm  is  order 
optimal.  It  solves  the  thinnest  path  problem  in  both  the  1- 
D  and  1.5-D  networks.  In  particular,  no  algorithm,  even 
with  complete  location  information  of  the  eavesdroppers, 
can  obtain  a  thinner  path  than  the  NBI  algorithm  which 
does  not  require  knowledge  of  eavesdropper  locations. 

In  a  broader  context,  our  techniques  of  using  directed 
crosses  and  exposed  disk  hypergraphs  in  establishing 
the  NP-completeness  of  the  problem  may  spark  new 
tools  for  complexity  studies  in  geometrical  hypergraphs 
and  graphs.  The  bounding  techniques  and  the  use  of 
sphere  packing  results  in  analyzing  the  performance  of 
the  two  approximation  algorithms  may  also  find  other 
applications  in  algorithmic  analysis. 

C.  Related  Work 

There  is  a  large  body  of  literature  on  security  issues 
in  wireless  ad  hoc  networks  (see,  for  example,  [4],  [5]). 
However,  the  thinnest  path  problem  has  not  been  studied 
in  the  literature  except  in  [6].  Chechik  et  al.  studied  the 
thinnest  path  (referred  to  as  the  secluded  path  in  [6]) 
and  the  thinnest  Steiner  tree  in  graphs.  They  showed 
that  the  problem  in  a  general  graph  is  NP-complete  and 
strongly  inapproximable.  They  proposed  an  algorithm 
with  an  approximation  ratio  of  \/A  +  3  where  A  is 
the  maximum  degree  of  the  graph.  They  further  studied 
the  problem  in  several  special  graph  models  including 
graphs  with  bounded  degree,  hereditary  graphs,  and 
planar  graphs.  However,  their  study  focuses  on  the 
problem  in  topological  graphs,  whereas  we  focus  on 
hypergraphs  and  geometric  graphs.  The  results  obtained 
in  [6]  do  not  apply  to  special  hypergraphs  satisfying 
certain  geometric  properties  that  result  naturally  from 
the  communication  problem  studied  in  this  paper.  This 
paper  also  includes  several  new  results  on  the  thinnest 
path  problem  under  the  graph  model.  Specifically,  we 
establish  the  NP-completeness  of  the  problem  in  2-D 
disk  graphs  and  3-D  unit  disk  graphs.  The  results  in  [6] 
and  this  work  thus  complement  each  other  to  provide 
a  more  complete  picture  of  the  thinnest  path  problem 
under  different  (hyper)graph  models. 

In  the  general  context  of  algorithmic  studies  in  hy¬ 
pergraphs,  Ausiello  et  al.  [7]  tackled  the  problem  of 


finding  the  /i-optimal  hyperpath  where  //  is  a  general 
measure  on  hyperpaths  that  satisfies  a  certain  monotone 
property.  They  established  the  NP-completeness  of  this 
problem  for  general  measures.  The  thinnest  path  problem 
can  be  seen  as  an  ^-optimal  traversal  problem  with  the 
measure  //  given  by  the  number  of  vertices  covered 
by  the  path.  Since  this  is  a  special  measure,  their  NP- 
completeness  result  developed  under  general  measures 
does  not  apply.  Furthermore,  in  many  applications,  the 
resulting  hypergraphs  have  certain  topological  and/or  ge¬ 
ometrical  properties,  and  the  computational  complexities 
under  these  special  models  require  separate  analysis. 

Another  realated  problem  is  the  shortest  path  problem 
in  hypergraphs,  which  remains  a  polynomial-time  prob¬ 
lem  as  its  counterpart  under  the  graph  model.  The  static 
shortest  hyperpath  problem  was  considered  by  Knuth  [8] 
and  Gallo  et  al.  [3]  in  which  Dijkstra’s  algorithm  for 
graph  was  extended  to  obtain  the  shortest  hyperpaths. 
Ausiello  et  al.  proposed  a  dynamic  shortest  hyperpath 
algorithm  for  directed  hypergraphs,  considering  only  the 
incremental  problem  (i.e.,  network  changes  contain  only 
edge  insertion  and  weight  decrease)  with  the  weights  of 
all  hyperedges  limited  to  a  finite  set  of  numbers  [9]. 
In  [10],  Gao  et  al.  developed  the  first  fully  dynamic 
shortest  path  algorithms  for  general  hypergraphs.  As 
discussed  earlier,  the  thinnest  path  problem  is  fundamen¬ 
tally  different  and  significantly  more  complex  than  the 
shortest  path  problem. 

The  widest  path  problem  has  been  well  studied  under 
the  graph  model  [11],  [12],  and  the  existing  results 
can  be  easily  extended  to  hypergraphs.  The  widest  path 
problem  asks  for  a  path  whose  minimum  edge  weight 
along  the  path  is  maximized.  In  other  words,  the  width 
of  a  path  is  given  by  the  minimum  edge  weight  on 
that  path,  which  is  different  from  the  definition  of  path 
width  in  the  thinnest  path  problem  studied  in  this  paper. 
As  a  consequence,  the  widest  path  problem  is  not  the 
complement  of  the  thinnest  path  problem.  Since  the  tree 
structure  is  preserved  in  the  widest  path  problem  (i.e., 
the  widest  path  to  a  node  must  go  through  the  widest 
path  to  one  of  its  neighbors),  it  remains  a  polynomial 
time  problem.  The  thinnest  path  problem,  however,  is 
NP-complete  in  general. 

II.  Problem  Formulation 
A.  Basic  Concepts  of  Directed  Hypergraphs 

A  directed  hypergraph  H  =  ( V ,  E )  consists  of  a  set  V 
of  vertices  and  a  set  E  of  directed  hyperedges  [3]1.  Each 
directed  hyperedge  e  £  E  has  a  single  source  vertex  se 
and  a  non-empty  set  of  destination  vertices  Te.  We  let 
n  =  \V\  denote  the  number  of  vertices. 

A  disk  hypergraph  is  a  special  directed  hypergraph 
whose  topology  is  determined  by  a  set  of  points  V  = 
{tii, . . . ,  vn}  located  in  a  d-dimensional  Euclidean  space 
and  a  maximum  range  Ri  associated  with  each  vertex 
Vi.  There  exists  a  hyperedge  e  from  source  se  =  vt  to 
destination  set  Te  if  and  only  if  Te  consists  of  vertices 
located  within  the  d-dimensional  sphere  centered  at  Vi 

[n  [3],  it  is  referred  to  as  the  forward  hyperarcs. 
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Fig.  1:  (a):  Directed  hypergraph;  (b):  Disk  hypergraph;  (c):  Unit  disk 
hypergraph. 


with  a  radius  r  £  (0,  Ri],  A  unit  disk  hypergraph 
(UDH)  is  a  disk  hypergraph  with  unit  maximum  range 
( Ri  =  1)  for  all  vertices.  Fig.  1  shows  examples  of  a 
directed  hypergraph,  a  disk  hypergraph,  and  a  unit  disk 
hypergraph. 

A  ring  hypergraph  is  a  generalized  disk  hypergraph 
where  associated  with  each  vertex  Vi  is  a  minimum  range 
ri  as  well  as  a  maximum  range  R, .  Hyperedges  rooted 
at  Vi  are  formed  by  spheres  centered  at  Vi  with  radii 
satisfying  rj  <  r  <  Ri.  It  is  easy  to  see  that  a  disk 
hypergraph  is  a  ring  hypergraph  with  ry  =  0,  a  disk 
graph  is  a  ring  hypergraph  with  r,  =  R,  for  all  i,  and 
a  unit  disk  graph  (UDG)  is  a  ring  hypergraph  with  r.L  = 
Ri  =  1  for  all  i. 


B.  The  Thinnest  Path  Problem 

Consider  a  wireless  ad-hoc  network  with  n  nodes 
located  in  a  (/-dimension  Euclidean  space.  Each  node 
can  choose  the  power,  within  a  maximum  value,  for  the 
transmission  of  each  message.  The  chosen  power,  along 
with  the  signal  propagation  model,  determines  the  set 
of  neighbors  that  can  hear  the  message.  The  maximum 
transmission  power  is  in  general  different  across  nodes. 
The  objective  is  to  find  a  path  between  a  given  source- 
destination  pair  that  involves  the  minimum  number  of 
nodes  hearing  the  message. 

As  discussed  in  Sec.  I- A,  we  formulate  the  problem 
using  a  directed  hypergraph.  Each  node  is  a  vertex. 
The  directed  hyperedges  rooted  at  a  node  are  given  by 
distinct  neighbor  sets  of  this  node  feasible  under  its 
maximum  transmission  power  and  the  signal  propagation 
model.  Under  a  general  nonisotropic  propagation  model, 
we  end  up  with  a  general  hypergraph.  The  only  property 
the  resulting  hypergraph  has  is  the  monotonicity  of 
the  hyperedge  set.  Specifically,  the  hyperedges  rooted 
at  each  node  can  be  ordered  in  such  a  way  (say, 
ei,e2,...,e;)  that  Tei  C  Tei+1  and  |Tei|  =  |Tei+1|  -  1. 
This  is  due  to  the  nature  of  wireless  broadcasting  where 
nodes  reachable  under  transmission  power  77  can  also 
be  reached  under  any  power  greater  than  77.  Under  an 
isotropic  propagation  model,  we  end  up  with  a  disk 
hypergraph.  If  all  nodes  have  the  same  maximum  range2, 
we  have  a  unit  disk  hypergraph.  This  hypergraph  model 
also  applies  to  networks  with  eavesdroppers.  Each  eaves¬ 
dropper  can  be  seen  as  a  node  with  zero  transmission 
range.  It  is  thus  a  vertex  with  no  outgoing  hyperedges. 

Given  a  source-destination  pair  (s,  t),  a  hyperpath 
from  s  to  t  is  defined  as  a  sequence  of  hyperedges 


^Transmission  range  and  transmission  power  are  used  interchange¬ 
ably. 


L  =  {ei, . . . ,  emj  such  that  sei  £  Tei_1  for  1  <  i  <m, 
sei  =  s,  and  t  £  TSm .  Define  the  cover  L  of  L  to  be  the 
set  of  vertices  in  L,  i.e., 

l=  ur=1Teiu{a 

The  width  W (L)  is  then  given  by 
W{L)t\L\. 

The  thinnest  path  problem  asks  for  a  hyperpath  from 
s  to  i  with  the  minimum  width.  Note  that  choosing 
a  hyperedge  e  =  {se,Te}  simultaneously  chooses  the 
relaying  node  se  and  its  transmission  power  (determined 
by  Te). 

III.  NP-Complete  Problems 

In  this  section,  we  show  that  the  thinnest  path  (TP) 
problem  is  NP-complete  in  several  special  geometric  hy¬ 
pergraphs  and  graphs.  This  implies  the  NP-completeness 
of  the  problem  in  general  directed  hypergraphs. 

A.  TP  in  2 -D  Disk  Hypergraphs 

In  this  subsection,  we  prove  the  NP-completeness 
of  the  thinnest  path  problem  in  2-D  disk  hypergraphs. 
While  a  stronger  result  is  shown  in  the  next  subsection, 
the  proof  of  this  result  provides  the  main  building  block 
for  the  proof  of  the  next  result. 

The  result  is  established  through  a  reduction  from 
the  maximum  dominating  set  (MDS)  [13]  problem.  The 
MDS  problem  asks  for  the  minimum  subset  of  vertices 
in  a  given  graph  such  that  every  vertex  in  the  graph  is 
either  in  the  subset  or  a  direct  neighbor  of  a  vertex  in  the 
subset.  The  following  theorem  formally  establishes  the 
polynomial  reduction  (denoted  by  <p)  from  MDS  to  TP 
in  2-D  disk  hypergraphs.  Since  the  thinnest  path  problem 
is  clearly  in  the  NP  space,  this  theorem  establishes  the 
NP-completeness  of  TP  in  2-D  disk  hypergraphs. 

Theorem  1:  MDS  <p  TP  in  2-D  disk  hypergrpahs. 

To  prove  Theorem  1,  consider  an  MDS  problem  in  an 
arbitrary  graph  G.  We  first  construct  a  general  directed 
hypergraph  Hi  based  on  G  such  that  a  thinnest  path 
in  Hi  leads  to  an  MDS  in  G.  The  main  challenge  in 
the  proof  is  to  show  that  Hi  is  realizable  under  a  2-D 
disk  model.  There  are  two  main  difficulties.  First,  line 
crossing  is  inevitable  when  we  draw  Hi  on  a  2-D  plane. 
The  implementation  of  hyperedges  that  cross  each  other 
needs  special  care  to  avoid  unwanted  overhearing  that 
may  render  the  reduction  invalid.  Second,  the  geometric 
structure  of  2-D  disk  hypergraphs  dictates  that  there  are 
at  most  5  vertices  (even  with  arbitrary  ranges)  that  can 
reach  a  common  sixth  vertex  but  not  each  other.  It  is  thus 
challenging  to  implement  a  vertex  with  up  to  n  incoming 
hyperedges  in  Hi  while  preserving  the  reduction. 

Our  main  approach  to  overcoming  the  above  diffi¬ 
culties  is  to  allow  directed  overhearing.  Specifically, 
messages  transmitted  along  one  hyperedge  may  be  heard 
by  vertices  implementing  another  hyperedge  in  Hi,  but 
not  vise  verse.  By  carefully  choosing  the  directions  of  the 
introduced  overhearing,  we  ensure  that  the  resulting  2-D 
disk  hypergraph  H2,  while  having  a  different  topological 
structure  from  Hi,  preserves  the  reduction  from  MDS 
in  G. 
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Another  challenge  in  constructing  H2  is  to  ensure 
the  polynomial  nature  of  the  reduction.  The  number 
of  additional  vertices  added  in  H2  needs  to  be  in  a 
polynomial  order  in  the  size  of  G.  This  often  limits  the 
use  of  reduced  transmission  ranges  as  a  way  to  avoid 
unwanted  overhearing:  exponentially  small  transmission 
ranges  may  require  exponentially  many  vertices  to  con¬ 
nect  two  fixed  points. 

A  detailed  proof  is  given  in  Appendix  A. 

B.  TP  in  2 -D  Unit  Disk  Hypergraphs 

We  now  establish  the  NP-completeness  of  TP  in  2-D 
unit  disk  hypergraphs  (UDH).  The  proof  builds  upon 
the  proof  of  Theorem  1.  The  only  difference  is  that 
when  implementing  the  general  directed  hypergraph  Hi, 
we  no  longer  have  the  freedom  of  choosing  the  maxi¬ 
mum  transmission  range  of  each  vertex.  This  presents 
a  non-trivial  challenge.  As  stated  in  Sec.  III-A,  our 
approach  to  circumventing  the  constraints  imposed  by 
the  geometrical  structures  of  2-D  disk  hypergraphs  is  to 
allow  directed  overhearing,  which  is  achieved  by  care¬ 
fully  choosing  different  maximum  transmission  ranges 
of  various  vertices.  To  implement  a  2-D  UDH  for  the 
reduction,  however,  all  vertices  must  have  the  same 
maximum  transmission  range. 

To  address  this  issue,  we  introduce  a  special  type  of 
disk  hypergraphs,  called  exposed  disk  hypergraplis,  and 
show  that  TP  in  k- D  exposed  disk  hypergraphs  can  be 
reduced  to  TP  in  k- D  UDH  for  any  k  >  2.  We  then  show 
that  the  2-D  disk  hypergraph  H2  constructed  in  the  proof 
of  Theorem  1  can  be  modified  to  an  exposed  hypergraph 
while  preserving  the  reduction.  We  thus  arrive  at  the  NP- 
completeness  of  TP  in  2-D  UDH  based  on  the  transitivity 
of  polynomial  time  reduction. 

Definition  1:  In  a  disk  hypergraph  H  =  (V,E),  let 
tv  denote  the  closest  non-neighbor  3  of  v.  Define4 

€v — 2  {d(,V,Tv')  Rv') : 

where  d(v,  tv)  is  the  distance  between  v  and  r„  (ev  is  set 
to  1  when  v  does  not  have  non-neighbors).  An  exposed 
area  <1>„  of  v  is  defined  as 

$v=Dv,Rv+ev\  |^J  DUtpu, 

UtzV 

where  Dv /r  denotes  the  closed  ball  centered  at  v  with 
radius  r.  A  disk  hypergraph  is  exposed  if  every  vertex 
has  a  non-empty  exposed  area. 


Fig.  2:  Exposed  hypergraphs  and  exposed  areas  (Hi  is  not  exposed 
since  v  has  an  empty  exposed  area;  H2  and  H3  are  exposed). 

Lemma  1:  TP  in  fc-D  exposed  disk  hypergraphs  <p 
TP  in  k-D  UDH. 

3  A  vertex  is  a  non-neighbor  of  v  if  it  is  outside  the  maximum  range 
Rv  of  v. 

4The  parameter  can  be  change  to  an  arbitrary  positive  value 
smaller  than  1. 


Proof:  The  basic  idea  is  to  place  super  vertices  at 
specific  locations  in  exposed  areas  to  force  vertices  on 
a  thinnest  path  to  use  transmission  ranges  smaller  than 
the  maximum  value.  The  problem  is  thus  transformed 
to  the  case  with  disk  hypergraphs  where  vertices  may 
have  different  maximum  transmission  ranges.  A  detailed 
proof  is  given  in  Appendix  B.  ■ 

With  Lemma  1  providing  a  bridge  between  disk  and 
unit  disk  hypergraphs,  all  we  need  to  show  is  that  MDS 
can  be  reduced  to  TP  in  2-D  exposed  disk  hypergraphs. 

Lemma  2:  MDS  <p  TP  in  2-D  exposed  disk  hyper- 
grpahs. 

Proof:  See  Appendix  C.  ■ 

Based  on  Lemma  1  and  Lemma  2,  we  arrive  at  the 
following  theorem. 

Theorem  2:  MDS  <P  TP  in  2-D  UDH. 

C.  TP  in  2-D  Disk  Graphs  and  3 -D  Unit  Disk  Graphs 
In  this  subsection,  we  consider  the  thinnest  path 
problem  in  disk  graphs  and  unit  disk  graphs  (UDG). 
Recall  that  disk  and  unit  disk  graphs  are  special  ring 
hypergraphs  with  r,  =  f?*  and  r,  =  R.,  =  1,  respectively. 
In  other  words,  they  can  be  seen  as  hypergraphs  where 
each  vertex  has  only  one  outgoing  hyperedge  directed  to 
its  prefixed  neighbor  set  (determined  by  its  fixed  trans¬ 
mission  power).  This  also  shows  that  disk  hypergraphs 
and  disk  graphs  are  not  special  cases  of  each  other.  Given 
the  same  set  of  vertices  and  their  associated  maximum 
ranges,  a  disk  hypergraph  has  a  topology  different  from 
that  of  a  disk  graph:  each  vertex  in  general  has  more 
than  one  outgoing  hyperedge  due  to  the  freedom  of 
using  smaller  transmission  ranges.  The  same  holds  for 
UDH  and  UDG.  As  a  consequence,  the  complexity  of 
TP  in  disk  and  unit  disk  graphs  cannot  be  inferred  from 
Theorems  1  and  2,  and  needs  to  be  studied  separately. 
Theorem  3:  MDS  <p  TP  in  2-D  disk  graphs. 

Proof:  In  the  proof  of  Theorem  1 ,  the  vertices  along 
the  thinnest  path  in  the  constructed  2-D  disk  hypergraph 
H2  all  use  their  maximum  ranges.  Thus,  MDS  in  G  can 
be  reduced  to  TP  in  a  disk  graph  constructed  from  H2 
by  including  only  those  hyperedges  associated  with  the 
maximum  range  of  each  vertex.  ■ 

Next  we  consider  TP  in  UDG.  Unfortunately,  the 
approach  through  exposed  disk  hypergraphs  used  in 
showing  the  NP-completeness  of  TP  in  UDH  does  not 
apply,  since  it  hinges  on  vertices  being  able  to  use  any 
transmission  range  smaller  than  a  maximum  value.  The 
difficulty,  however,  can  be  circumvented  for  3-D  UDG 
as  shown  in  the  following  theorem. 

Theorem  4:  MDS  in  degree-3  graphs  <p  TP  in  3-D 
UDG. 

The  proof  is  similar  to  that  of  Theorem  1  with  two 
main  differences.  First,  line  crosses  are  implemented  by 
using  the  third  dimension  to  “go  around”,  rather  than 
using  different  transmission  ranges  (a  luxury  absent  in 
UDG)  to  create  directed  crosses.  Second,  reduction  from 
MDS  in  graphs  with  a  maximum  degree  of  3  ensures 
that  there  are  at  most  4  incoming  edges  to  each  super 
vertex  in  the  reduced  UDG.  This  makes  the  geometric 
constraint  on  the  number  (at  most  11  in  a  3-D  Euclidean 
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space)  of  vertices  that  can  reach  a  common  vertex  but  not 
each  other  inconsequential5.  A  detailed  proof  is  given  in 
Appendix  D. 

Note  that  using  a  reduction  from  MDS  in  graphs  with 
a  constant  maximum  degree  rather  than  MDS  in  general 
graphs  leads  to  a  weaker  statement.  While  MDS  in  both 
cases  are  NP-complete,  the  former  is  approximable  with 
a  constant  ratio,  and  the  latter  a  ratio  of  O(logn). 
Theorems  1-3  thus  give  a  logn  order  lower  bound 
on  the  approximation  ratio  of  those  problems  whereas 
Theorem  4  provides  a  constant  lower  bound. 

IV.  Polynomial  Complexity  Problems 

In  this  section,  we  consider  the  thinnest  path  prob¬ 
lem  in  1-D  networks.  We  show  that  the  problem  is 
polynomial  time  by  constructing  an  algorithm  with  time 
complexity  of  O(n).  Since  the  input  data  has  size 
O  (n  'j ,  the  proposed  algorithm  is  order-optimal.  We  then 
consider  the  1.5-D  problem  and  show  that  the  algorithm 
developed  for  1-D  networks  directly  applies  to  the  1.5-D 
problem. 

A.  1  -D  Networks 

Consider  a  network  under  a  general  propagation  model 
with  n  nodes  located  on  a  straight  line.  Each  vertex 
Vi  is  associated  with  a  coordinate  xt  on  the  line  (the 
vertex  index  Vi  and  its  location  x%  are  often  used 
interchangeably).  Without  loss  of  generality,  we  assume 
that  x\  <  X2  <  . . .  <  xn. 


Fig.  3:  A  1-D  network  (circles  represent  maximum  ranges  under  a 
disk  propagation  model). 


It  is  clear  that  every  node  located  between  the  source 
s  and  the  destination  t  (see  Fig  3)  will  hear  the  message 
no  matter  which  path  is  chosen  and  all  nodes  to  the  right 
of  t  can  be  excluded  from  the  thinnest  path.  Therefore, 
finding  the  thinnest  path  is  equivalent  to  minimizing  the 
number  of  vertices  to  the  left  of  s  that  can  overhear 
the  message.  The  problem  is  nontrivial.  Due  to  the 
arbitrariness  of  the  node  locations  and  propagation  range, 
a  forward  path  (i.e.,  every  hop  moves  the  message  to  the 
right  toward  t)  from  s  to  t  may  not  exist  and  nodes  to 
the  left  of  s  may  need  to  act  as  relays.  The  question  is 
thus  how  to  efficiently  find  out  whether  a  forward  path 
exists  and  if  not,  which  set  of  nodes  to  the  left  of  s  need 
to  relay  the  message. 

We  propose  an  algorithm  based  on  nested  backward 
induction  (NBI).  For  each  vertex  v,  we  define  its  prede¬ 
cessor  pv  to  be  the  nearest  vertex  on  the  left  side  of  v 

5We  can  consider  a  reduction  from  MDS  in  graphs  with  a  maximum 
degree  up  to  9  (see  Appendix  D). 


that  can  reach  v: 

pv  =argmax{x„  :  xu  <  xv, 

u£V 

3 e  €  E  s.t.  se  =  u  and  v  £  Te}.  (1) 
Thus,  in  order  to  reach  v,  its  predecessor  pv  or  a  vertex 
to  the  left  of  pv  has  to  transmit.  In  other  words,  those 
vertices  between  pv  and  v  cannot  directly  reach  v. 
Equivalently,  any  vertex  to  the  right  of  v  can  only  hear 
a  message  from  s  through  a  relay  by  pv  or  a  vertex  to 
the  left  of  pv. 

The  NBI  algorithm  is  then  carried  out  in  two  steps. 
In  the  first  step,  the  predecessors  of  certain  vertices  are 
obtained  one  by  one  starting  from  t  moving  toward  s. 
Specifically,  the  predecessor  of  t,  denoted  by  ui  =  pt,  is 
first  obtained.  If  xUl  <  xs,  then  the  first  step  terminates. 
Otherwise,  the  predecessor  of  u\,  denoted  by  U2  =  pUl, 
is  obtained  and  its  location  compared  with  xs.  The 
same  procedure  continues  until  the  currently  obtained 
predecessor  is  to  the  left  of  s  or  is  s  itself.  The  first  step 
thus  produces  a  sequence  of  vertices  m,  U2,  ■  ■  ■ ,  ui  with 
ui  =  pt ,  u2  =  pUl,  •  •  ■,  ui  =  pui_1  and  xui  <  xs.  Then 
L\  :=  . . .  ,ui,t}  is  a  valid  path  from  ui  to 

t.  If  ui  =  s,  the  algorithm  terminates,  and  the  thinnest 
path  from  s  to  t  is  given  by  L  \ .  Otherwise,  we  carry 
out  Step  2  of  the  algorithm  where  we  find  a  path  from 
s  to  ui.  Specifically,  let  V'  denote  the  set  of  vertices 
located  between  ui  and  ui-±  including  ui  but  not  ui-\. 
Fet  E'  denote  the  set  of  all  hyperedges  whose  source  and 
destination  vertices  are  in  V' .  As  shown  in  Appendix  E 
on  the  correctness  of  the  algorithm,  any  hyperpath  L2 
from  s  to  ui  in  the  sub-hypergraph  H'  =  (V',E’) 
concatenated  with  L\  gives  a  thinnest  path  from  s  to 
t.  Finding  such  an  L2  can  be  easily  done  by  a  breadth- 
first  search  (BFS)  in  H' .  However,  the  resulting  time 
complexity  is  0(n2).  Hence,  we  propose  a  special  BFS 
procedure  that  reduces  the  time  complexity  to  O  (n). 
The  trick  here  is  to  set  up  two  pointers,  ki  and  kr ,  to 
the  locations  of  the  leftmost  and  the  rightmost  vertices 
in  V'  that  have  been  discovered.  Due  to  the  geometric 
structure  of  the  1-D  network,  each  time  we  only  need  to 
search  vertices  to  the  left  of  ki  and  vertices  to  the  right 
of  kr.  The  detailed  algorithm  is  given  below. 

1.  Enqueue  s,  set  ki  and  kr  to  the  index  of  s. 

2.  Repeat  until  the  queue  is  empty  or  ui  is  found: 

-  Dequeue  a  vertex  v  and  examine  it 

-  If  v  =  ui,  go  to  step  4. 

-  Otherwise, 

*  While  v  can  reach  Vkr+i 

■  Enqueue  Vkr+ 1  and  kr  4—  kr  +  1 

•  Set  the  parent  of  Vkr+ 1  to  v 

*  While  v  can  reach  v^-i 

•  Enqueue  Vk,~  1  and  fc;  fcz  —  1 

•  Set  the  parent  of  v^-i  to  v 

3.  If  the  Queue  is  empty,  return  “no  path  from  s  to  t”. 

4.  Trace  back  to  s  and  return  L2- 

The  following  theorem  establishes  the  correctness  of 
the  proposed  NBI  algorithm.  Furthermore,  it  reveals 
a  strong  property  of  the  path  obtained  by  NBI  under 
a  disk  propagation  model.  Specifically,  under  a  disk 
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propagation  model,  we  define  the  covered  area  A{L) 
of  a  hyperpath  L  =  {ei, . . . ,  emj  as 

m 

A(L)=  (J  DSei!rH,  (2) 

i= 1 

where  rei  is  the  minimum  transmission  range  that  in¬ 
duces  hyperedge  e*,  i.e., 

r€i  =  ma x{d(sei,i>)}.  (3) 

Theorem  5  shows  that  the  covered  area  of  the  path 
obtained  by  NBI  is  a  subset  of  the  covered  area  of  any 
feasible  path  from  s  to  t. 

Theorem  5:  NBI  algorithm  finds  the  thinnest  path 
L* .  Furthermore,  under  a  disk  propagation  model,  given 
any  valid  path  L  from  s  to  t,  we  have  Al*  C  A^. 

Proof:  See  Appendix  E.  ■ 

Theorem  6:  The  time  complexity  of  the  NBI  algo¬ 

rithm  is  O  (n) . 

Proof:  The  O  (n)  complexity  of  the  first  step  of  NBI 
is  readily  seen.  In  the  second  step,  the  time  complexity 
is  dominated  by  updating  the  queue  at  each  iteration.  Let 
k  denote  the  number  of  iterations  in  step  2.  Note  that  we 
only  check  m*  +  2  vertices  at  iteration  i,  where  m,  is 
the  number  of  new  vertices  that  have  been  enqueued  at 
this  iteration  and  l  —  \V\-  Also  k  is  bounded 
by  \V'\.  Hence  the  total  time  complexity  of  this  step  is 
bounded  by  +  2)  <  3|H'|.  We  thus  arrive  at 

the  theorem.  ■ 

B.  1.5 -D  Networks 

We  now  consider  the  1.5-D  problem  where  in-network 
nodes  are  located  on  a  line  and  eavesdroppers  are  located 
in  a  ('/-dimensional  space  that  contains  the  line  network. 
We  focus  on  the  disk  propagation  model.  A  unit  cost 
is  incurred  for  each  in-network  node  that  hears  the 
message  and  a  non-negative  cost  c  is  incurred  for  each 
eavesdropper  that  hears  the  message.  The  objective  is  to 
find  a  path  L*  from  s  to  t  with  the  minimum  total  cost: 

L*=  arg  min  {  V''  c(v)}  (4) 

L={ei  .,em}  ' 

1  1  veA(L) 

where  c( v)  is  the  cost  for  vertex  v,  and  A(L)  is  the 
covered  area  of  path  L  as  defined  in  (2). 

Based  on  Theorem  5,  it  is  easy  to  see  that  NBI 
provides  the  optimal  solution  to  the  1.5-D  thinnest 
path  problem  without  knowledge  of  the  eavesdroppers 
locations.  More  specifically,  no  algorithm,  even  with 
complete  knowledge  of  the  locations  of  the  eavesdrop¬ 
pers,  can  obtain  a  thinner  path  than  NBI  which  does  not 
require  location  knowledge  of  the  eavesdroppers. 

V.  Approximation  Algorithms 

In  this  section,  we  introduce  two  approximation  algo¬ 
rithms  for  the  thinnest  path  problem  and  analyse  their 
performance  in  different  types  of  hypergraphs. 

A.  Shortest  Path  Based  Approximation  Algorithm 

Given  a  general  directed  hypergraph  H  with  source 
vertex  s  and  destination  vertex  t,  we  set  the  weight  of 
a  hyperedge  to  be  the  number  of  destination  vertices  in 
this  hyperedge: 

w(e)=\Te\  (5) 


The  shortest  hyperpath  algorithm  from  s  to  t  is  then  ob¬ 
tained  under  this  weight  definition  as  an  approximation 
of  the  thinnest  path.  The  following  theorem  quantifies 
the  performance  of  this  shortest  path  based  algorithm 
(SPBA). 

Theorem  7:  The  SPBA  algorithm  provides  a  \/j- 
approximation  for  TP  in  general  directed  hypergraphs, 
a  2(1  +  2a)  -approximation  for  d— dimensional  ring 
hypergraphs  with  a  =  — r — max,'igv  R' -jj — yt- 

°  r  max{minVigy  r*j,  minUiUgy  a{u,v)\ 

Additionally,  the  ratio  ^/lj  of  the  SPBA  algorithms  is 
asymptotically  tight  even  in  2-D  disk  hypergraphs. 

Proof:  See  Appendix  F.  ■ 

B.  Tree  Structure  Based  Approximation  Algorithm 

Approximation  occurs  in  two  places  in  SPBA.  First, 
the  width  of  a  path  is  approximated  by  the  sum  of 
the  widths  of  the  hyperedges  on  that  path.  Second,  the 
thinnest  path  to  a  vertex  is  assumed  to  go  through  the 
thinnest  path  to  one  of  its  incoming  neighbors.  The  first 
approximation  can  be  avoided  while  maintaining  the 
polynomial  nature  of  the  approximation  algorithm.  In 
particular,  we  can  ensure  that  the  width  of  a  path  is  cor¬ 
rectly  obtained  by  using  the  set  union  operation  instead 
of  summation.  The  assumption  on  the  tree  structure  of 
the  thinnest  paths  allows  us  to  use  Dijkstra’s  algorithm 
with  some  modifications.  Specifically,  for  each  vertex, 
we  need  to  store  the  current  thinnest  path  from  s  to  this 
vertex  rather  than  only  the  width  of  this  path  and  the 
parent  of  this  vertex  on  this  path.  This  allows  us  to  take 
the  set  union  operation  when  we  update  the  neighbors 
of  this  vertex.  Given  below  is  the  performance  of  this 
tree  structure  based  algorithm  (TSBA). 

Theorem  8:  The  TSBA  algorithm  provides  a  2f_1- 
approximation  for  general  directed  hypergraphs,  2(1  + 
2a)d-approximation  for  d— dimensional  ring  hyper¬ 
graphs  with  a  =  j — t — max«»gi  Ri — yt-  Addi- 

°  ^  max{min^.gv  aiu,!/)} 

tionally,  the  ratio  -  of  the  TSBA  algorithm  is  tight 
in  general  directed  hypergraphs  and  asymptotically  tight 
in  disk  hypergraphs  in  the  worst  case. 

Proof:  See  Appendix  G.  ■ 


C.  Performance  Comparison 

While  the  approximation  ratio  of  TSBA  is  better  than 
that  of  SPBA,  these  are  worst-case  performances  and  do 
not  imply  that  TSBA  outperforms  SPBA  in  every  case 
as  shown  in  Fig  4. 


Fig.  4:  An  example  where  SPBA  outperforms  TSBA  (SPBA  returns 
the  path  that  goes  through  all  solid  black  hyperedges  to  v  and  then  to 
t,  which  is  the  thinnest  path;  TSBA  returns  the  path  that  contains  all 
dashed  hyperedges,  which  is  not  the  thinnest  path). 

Fig.  5  shows  the  average  performance  of  these  two 
algorithms.  We  see  that  both  algorithms  have  relatively 
small  approximation  ratios  growing  sub-linearly  with  the 
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Fig.  5:  Average  performance  of  SPBA  and  TSBA  (a  2-D  network  with 
n  vertices  uniformly  and  randomly  distributed  on  a  ^  x  ^  square  with 
p  =  1.5;  the  maximum  range  of  each  vertex  randomly  chosen  from 
interval  [Rmin,  Rmax]  with  Rmin  =  l,  Rmax  =  5;  average  taken 
over  1000  such  random  2-D  disk  hypergraphs). 

number  of  vertices.  In  general,  TSBA  outperforms  SPBA 
on  average,  as  also  demonstrated  in  a  number  of  other 
simulation  results  (omitted  due  to  the  space  limit). 

VI.  Conclusion 

We  studied  the  complexity  and  developed  optimal  and 
approximation  algorithms  for  the  thinnest  path  problem 
for  secure  communications  in  wireless  ad  hoc  networks. 
In  establishing  the  NP-completeness  of  the  problem,  our 
techniques  of  using  directed  crosses  and  exposed  disk 
hypergraphs  may  spark  new  tools  for  complexity  studies 
in  geometrical  hypergraphs  and  graphs.  The  bounding 
techniques  and  the  use  of  sphere  packing  results  in 
analyzing  the  performance  of  the  two  approximation  al¬ 
gorithms  may  also  find  other  applications  in  algorithmic 
analysis. 

Appendix  A.  Proof  of  Theorem  1 

A.  Reduction  from  MDS  to  TP  in  A  General  Directed 
Hypergraph  1 1\ 

Consider  the  MDS  problem  in  a  graph  G  with  n 
vertices  Vi, ...  ,vn.  We  construct  a  directed  hypergraph 
H i  based  on  G  as  follows.  The  vertex  set  of  H\  includes 
the  n  vertices  of  G  augmented  by  a  destination  vertex 
vn+\  and  n  super  vertices  vf , . . . ,  u* .  A  super  vertex  v? 
corresponds  to  the  normal  vertex  v,  and  is  a  set  of  ns 
normal  vertices.  The  hyperedges  in  Hi  are  all  rooted  at 
the  normal  vertices  v\, ...  ,vn.  Specifically,  rooted  at  v.t 
(1  <  i  <  n)  are  fc,  +  1  directed  hyperedges,  where  ki 
is  the  degree  of  Vi  in  G.  Each  hyperedge  rooted  at  v.t 
has  two  destinations:  Vi+i  and  a  super  vertex  vf  whose 
corresponding  normal  vertex  Vj  dominates6  v,  in  the 
original  graph  G.  Fig.  6  is  an  example  illustrating  the 
construction  of  Hi  from  G. 

From  the  construction  of  Hi,  we  see  that  any  path 
from  vi  to  vn+i  must  traverse  through  all  normal  vertices 
one  by  one.  There  are  multiple  hyperedges  leading 
from  Vi  to  Uj+i,  each  involving  a  super  vertex  that 
corresponds  to  a  dominating  node  of  Vi  in  G.  Thus, 
choosing  a  hyperedge  going  from  Vi  to  vl+\  is  equivalent 
to  choosing  a  dominating  node  of  Vi  in  G.  Since  every 
path  from  Vi  to  vn+\  includes  all  the  n  +  1  normal 
vertices,  the  thinnest  path  is  given  by  the  one  with  the 

6A  vertex  in  a  graph  is  dominated  by  itself  and  any  of  its  one-hop 
neighbors. 


minimum  number  of  super  vertices,  thus  leading  to  the 
MDS  in  G.  At  this  point,  the  size  ns  of  a  super  vertex 
can  be  any  positive  integer.  As  will  become  clear  later, 
to  implement  Hi  under  a  2-D  disk  model,  additional 
normal  vertices  need  to  be  added.  As  a  consequence, 
paths  from  iq  to  vn+i  may  include  different  numbers  of 
normal  vertices.  To  preserve  the  reduction,  we  need  to 
make  sure  that  the  width  of  a  path  is  dominated  by  the 
number  of  super  vertices  it  covers.  This  can  be  achieved 
by  choosing  an  ns  sufficiently  large  (see  Appendix  A-D) 


Vi  V2  v3  V4  V5  V6 


Fig.  6:  The  construction  of  H i  from  G\  (a)  the  graph  G ;  (b)  the 
hypergraph  H\  (?;i  is  dominated  by  v\  and  vs  in  G.  We  thus  have 
two  hyperedges  rooted  at  v\  in  H±:  one  reaches  (v2,vf),  the  other 
(V2,V%).). 

The  following  lemma  formally  establishes  the  correct¬ 
ness  of  the  reduction. 

Lemma  3:  There  is  a  dominating  set  with  size  k  in 
G  if  and  only  if  there  is  a  path  from  vi  to  vn+  i  in  Hi 
with  width  kns  +  n  +  1. 

Proof:  First,  assume  that  G  has  a  dominating  set  S 
with  size  k.  By  the  definition  of  dominating  set,  for  each 
vertex  tq  in  G,  there  is  a  vertex  Vj  €  S  that  dominates  rq. 
From  the  construction  of  Hi,  there  exists  a  hyperedge  e, 
(i  =  1, . . . ,  n)  in  JTi  directed  from  tq  to  vertex  Vi+i  and 
super  vertex  v?  corresponding  to  the  dominating  node  Vj 
in  G.  Thus,  the  hyperpath  {ei, . . .  ,enj  is  a  path  from 
Vi  to  vn+i  with  width  kns  +  n  +  1.  The  width  comes 
from  the  fact  that  all  n  +  1  vertices  in  Vh1  are  on  the 
path  along  with  k  super  vertices,  each  consisting  of  ns 
normal  vertices. 

Conversely,  assume  that  there  exists  a  path  from  v\ 
to  vn+i  in  Hi  with  width  kns  +  n  +  1.  Based  on  the 
construction  of  H,  every  path  from  vi  to  vn+\  consists 
of  n  hyperedges  rooted  at  each  of  the  n  normal  vertices 
Vi, ...  ,vn.  Thus,  a  path  with  width  kns  +  n  +  1  must 
contain  k  super  vertices.  From  the  construction  of  the 
hyperedges,  we  conclude  that  the  vertices  in  G  that 
correspond  to  those  k  super  vertices  along  the  given  path 
form  a  dominating  set  with  size  k.  ■ 

B.  A  2-D  Grid  Representation  of  Hi 

The  directed  hypergraph  Hi  obtained  above  does  not 
satisfy  the  geometric  properties  of  2-D  disk  hypergraphs 
(see  Sec.  II).  To  prove  Theorem  1,  we  need  to  modify 
Hi  to  a  2-D  disk  hypergraph  H2  while  preserving 
the  reduction  from  MDS  in  G.  Our  approach  is  to 
realize  the  topological  structure  of  each  hyperedge  in 
Hi  by  adding  additional  vertices  with  carefully  chosen 
locations  and  maximum  ranges  to  lead  from  the  source 
vertex  to  the  destination  vertices  of  this  hyperedge. 
The  number  of  additional  vertices,  however,  should  be 
kept  at  a  polynomial  order  with  the  problem  size  to 
ensure  the  polynomial  nature  of  the  reduction.  This  can 
be  achieved  by  adding  vertices  on  a  2-D  grid  with  a 


constant  grid  spacing,  which  allows  a  constant  maximum 
range,  thus  polynomially  many  additional  vertices.  The 
detailed  implementation  of  Hi  under  a  2-D  disk  model 
is  given  in  the  next  subsection.  As  a  preparatory  step, 
we  show  in  this  subsection  that  the  hyperedges  in  Hi 
can  be  represented  by  line  segments  of  a  2-D  grid  with 
a  constant  grid  spacing. 

We  first  embed  the  normal  vertices  of  Hi  evenly 
in  a  horizontal  line  in  a  2-D  space  (see  Fig.  7  for 
an  illustration).  Below  this  line  is  a  2 n2  x  4 n2  unit 
grid.  There  are  4 n  vertical  lines  between  vt  and  vl+-\ 
(1  <  i  <  n)  that  are  partitioned  into  three  zones  (Cf , 
Cf,  Cf)  of  n,  2 n,  and  n  vertical  lines,  respectively.  The 
super  vertices  are  embedded  evenly  on  a  horizontal  line 
below  the  grid.  The  horizontal  position  of  super  vertex 
vf  is  between  Vi  and  Uj+i. 

Next,  we  specify  how  a  hyperedge  traverses  the  grid 
from  its  source  vertex  to  its  destination  vertices.  Recall 
that  every  hyperedge  in  Hi  is  directed  from  a  normal 
vertex  vt  to  a  super  vertex  and  the  next  normal 
vertex  u,+i.  To  preserve  the  reduction,  we  need  to  ensure 
that  each  hyperedge  can  only  reach  its  normal  vertex 
destination  after  reaching  its  super  vertex  destination.  To 
facilitate  the  implementation  around  the  super  vertices 
(see  Appendix  A-C2),  we  designate  the  middle  zone 
Cf  between  vt  and  Vi+i  for  traveling  down  to  super 
vertex  vf  and  then  up  to  the  corresponding  normal  vertex 
destination  (see  region  Cf  in  Fig.  7).  Each  hyperedge 
involving  vf  has  two  designated  vertical  lines  in  Cf 
(one  for  going  down  to,  the  other  going  up  from  the 
super  vertex).  To  connect  the  designated  vertical  lines 
in  zone  Cf  with  the  source  vertex  and  then  to  the 
normal  destination  vertex,  we  designate  two  horizontal 
lines  for  each  hyperedge.  The  traverse  of  the  hyperedge 
is  completed  by  designating  one  vertical  line  in  Cf 
and  one  in  Cf  to  connect  the  normal  vertices  with  the 
corresponding  designated  horizontal  lines.  Since  there 
are  at  most  n2  hyperedges,  the  designed  grid  size  is 
sufficient  to  ensure  that  each  hyperedge  traverses  through 
a  distinct  set  of  line  segments  in  the  grid.  . 


Fig.  7:  A  2-D  grid  representation  of  Hi  (the  two  hyperedges  rooted 
at  v  i  from  the  example  given  in  Fig.  6  are  illustrated  in  green  and 
blue,  respectively). 


C.  Implementing  Hi  under  A  2-D  Disk  Model 

Based  on  the  2-D  grid  representation  of  Hi,  we  can 
construct  a  2-D  disk  hypergraph  Hi  that  preserves  the 
reduction.  Specifically,  we  place  a  sequence  of  evenly 
spaced  normal  vertices  with  a  constant  maximum  range 


along  the  line  segments  in  the  grid  that  form  each 
hyperedge  of  H\.  The  distance  between  two  adjacent 
vertices  is  set  to  their  maximum  range.  The  constant 
maximum  range  can  be  set  sufficiently  small  (say,  A)  to 
avoid  overhearing  across  vertices  on  different  hyperedges 
that  may  render  the  reduction  invalid.  There  are  two 
issues  that  remain  to  be  addressed:  the  implementation 
of  crosses  and  that  around  super  vertices. 

1 )  Implementaion  of  Crosses:  The  line  crossing  in 
the  grid  representation  of  Hi  make  overhearing  across 
hyperedges  inevitable.  However,  by  exploiting  the  free¬ 
dom  of  choosing  the  maximum  range  for  each  vertex,  we 
can  implement  directed  crosses  that  allow  us  to  preserve 
the  reduction.  Specifically,  when  two  line  segments  in  the 
grid  representation  cross,  we  can  choose  the  maximum 
ranges  of  the  vertices  along  these  two  lines  in  such  a 
way  that  messages  transmitted  over  one  line  can  be  heard 
by  vertices  on  the  other  but  not  vise  verse.  A  specific 
implementation  is  given  in  Fig.  8. 


Fig.  8:  A  disk  hypergraph  implementation  of  a  directed  cross  where  the 
circles  represent  the  maximum  range  of  vertices  (messages  transmitted 
on  the  blue  line  can  be  heard  by  nodes  on  the  red  line,  but  not  vise 
versa). 

Next,  we  show  how  carefully  choosing  the  direction 
of  each  cross  allows  us  to  preserve  the  reduction.  The 
cross  directions  are  defined  by  assigning  a  level  index  to 
each  line  segment  in  the  grid  representation.  Specifically, 
For  a  hyperedge  rooted  at  Vi  in  Hi,  its  line  segments 
before  and  after  reaching  the  super  vertex  destination 
have  levels  i  and  i  + 1,  respectively.  Then,  each  cross  has 
a  direction  pointing  from  the  higher  level  segment  to  the 
lower  one  (i.e.,  messages  transmitted  on  the  higher  level 
segment  can  be  heard  by  the  vertices  along  the  lower 
level  segment  but  not  vise  versa).  If  the  two  segments 
have  the  save  level,  the  direction  of  the  cross  can  be 
arbitrary.  To  see  that  this  directed  implementation  of 
crosses  preserves  the  reduction,  we  only  need  to  notice 
that  any  path  from  Vi  to  still  needs  to  go  through 
all  the  n  normal  vertices  one  by  one  and  must  reach  a 
super  vertex  before  reaching  the  next  normal  vertex. 

2)  Implementation  Around  Super  Vertices:  Recall 
that  a  super  vertex  in  Hi  is  a  set  of  ns  normal  vertices 
that  have  no  outgoing  hyperedges.  It  can  be  implemented 
by  ns  points  with  zero  maximum  range  and  located 
sufficiently  close  to  each  other  (so  that  any  path  from 
Vi  to  vn+i  in  Hi  includes  either  all  of  them  or  none  of 
them). 

Consider  first  the  implementation  of  one  incoming 
hyperedge  to  a  super  vertex  vf .  Recall  that  in  the  2- 
D  grid  representation  of  Hi,  a  hyperedge  approaches 
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and  leaves  v}  through  two  vertical  lines  in  zone  C2 
(see  Fig.  7).  One  implementation  of  this  U-turn  around 
Vj  is  to  add  6  normal  vertices  with  specific  maximum 
ranges  and  locations.  As  shown  in  Fig.  9,  these  6  vertices 
include  three  anchor  vertices  u j~,  u and  uf  with 
maximum  range  r,  two  interface  vertices  vf  and  v3  that 
connect  with  the  grid,  and  a  bridging  vertex  /in,  all  with 
maximum  range  |.  The  value  of  r  and  the  connection 
with  the  grid  will  be  specified  later. 


Fig.  9:  Implementation  of  one  hyperedge  passing  through  a  super 
vertex.  Starting  from  uf .  the  message  traverse  to  uf  through  /in, 
u3,  ui>  uf.  The  super  vertex  hears  the  message  in  the  transmission 
from  nj  t°  tq  . 

A  challenge  remains  in  the  implementation  of  up  to  n 
incoming  hyperedges  to  the  same  super  vertex.  Note  that 
under  a  2-D  disk  model,  one  can  at  most  have  5  vertices 
(even  with  arbitrary  ranges)  that  reach  a  common  sixth 
vertex  but  not  each  other.  The  key  to  circumvent  this 
difficulty  is  to  allow  directed  overhearing,  similar  to  the 
idea  behind  the  implementation  of  the  crosses.  Specifi¬ 
cally,  the  reduction  is  preserved  as  long  as  a  hyperedge 
rooted  at  v3  cannot  overhear  a  message  transmitted  over 
a  hyperedge  rooted  at  for  any  i  <  j.  The  detailed 
implementation  is  as  follows.  The  fist  step  is  to  designate 
the  vertical  lines  in  zone  C2  to  the  incoming  hyperedges 
of  Vj  based  on  the  indices  of  their  source  vertices. 
Specifically,  the  incoming  hyperedge  with  the  smallest 
source  vertex  index  takes  the  two  center  most  lines  in 
Cj,  and  so  on.  Consider  first  the  implementation  of 
the  two  incoming  hyperedges  (say,  e\  and  e2)  with  the 
smallest  source  vertex  indices.  As  shown  in  Fig.  10,  we 
first  implement  e\  as  described  above  (see  Fig.  9).  The 
structure  of  the  implementation  of  e2  is  similar  except 
that  the  maximum  range  of  the  anchor  vertices  uf  ,  u 
and  ulf  is  set  to  4r  to  prevent  unwanted  overhearing. 
As  a  consequence,  more  bridging  vertices  (/i2i,  P22,  P23 
with  maximum  range  r,  and  2 r,  respectively)  are 
needed  to  connect  the  interface  vertex  u.J  to  the  anchor 
vertex  uf  .  Note  that  no  vertices  along  e2  (the  centers  of 
the  blue  circles  in  Fig.  10)  are  in  the  range  of  any  vertices 
along  e\  (the  green  circles).  The  correct  direction  of 
overhearing  is  thus  ensured. 

The  same  procedure  continues  for  any  additional 
incoming  hyperedges  to  vf,  in  the  ascending  order  of 
their  source  vertex  indices  in  H\.  Note  that  the  range  of 
the  anchor  vertices  in  the  fcth  hyperedge  is  4fer,  growing 
exponentially  with  k.  The  maximum  ranges  (specifically, 
5,  r, . . . ,  4 yr)  of  the  bridging  vertices  {fiki}  are  chosen 
to  preserve  the  polynomial  nature  of  the  reduction.  In  this 
way,  the  number  of  additional  vertices  for  implementing 
the  Arth  hyperedge  is  2k  +  4,  and  the  total  number  of 
additional  vertices  around  one  super  vertex  is  at  most 


n2  +  3  n. 


Fig.  10:  Implementation  of  the  second  incoming  hyperedge  to  a  super 
vertex. 


Next  we  consider  the  value  of  r  which  should  be  set 
sufficiently  small  to  avoid  overhearing  across  hyperedges 
leading  to  different  super  vertices.  Note  that  the  width 
of  the  area  covered  by  the  additional  vertices  around  a 
super  vertex  is  4  times  the  largest  maximum  range  of 
the  anchor  vertices.  We  thus  set  r  =  4 ~nn,  considering 
the  distance  between  two  adjacent  super  vertices  being 
4  n. 

The  last  issue  is  to  connect  the  interface  vertices  with 
the  grid.  Each  interface  vertex  needs  to  be  connected 
with  a  designated  vertical  line  in  C2.  While  the  vertical 
lines  in  C2  are  evenly  spaced,  the  horizontal  positions 
of  the  interface  vertices  have  an  exponential  structure 
due  to  the  exponentially  growing  range  of  the  anchor 
vertices.  Furthermore,  the  vertices  realizing  the  vertical 
lines  in  the  grid  have  a  constant  range,  whereas  the 
interface  vertices  have  an  exponentially  smaller  range 
of  r  =  4 ~nn.  If  we  connect  them  using  a  sequence 
of  vertices  with  a  constant  range,  unwanted  overhearing 
will  occur  near  the  interface  vertices.  On  the  other  hand, 
connecting  them  using  vertices  with  range  r  results 
in  an  exponential  number  of  additional  vertices.  To 
preserve  the  correctness  and  the  polynomial  nature  of  the 
reduction,  we  propose  the  scheme  detailed  in  Fig.  11. 

Since  the  generated  sequence  of  circles  {Ei}  are 
within  the  boundary  given  by  lines  AC  and  BD  and 
the  boundary  lines  corresponding  to  different  interface 
vertices  do  not  cross  (see  Fig.  12),  the  above  scheme 
does  not  introduce  overhearing,  thus  preserving  the  re¬ 
duction.  The  polynomial  nature  of  the  reduction  can  be 
shown  based  on  the  following  lemma. 

Lemma  4:  Consider  the  geometrical  scheme  de¬ 
scribed  in  Fig.  11.  Assume  ZCAB  >  j.  The  num¬ 
ber  of  circles  {Ei},  denoted  by  k,  satisfies  k  < 
2(1°Sfl11-fl2gfl2)L  +  1  when  ±  r2>  and  k  <  +  1 

when  R\  =  R2,  where  If ,  R2  denote  the  radii  of  circles 
0\  and  0‘2,  and  L  the  distance  between  lines  AB  and 
CD. 

Proof:  Assume  first  R\  f  If  -  Without  loss  of 
generality,  assume  If  >  f?2-  Since  line  AB  and  CD 
are  parallel,  the  three  lines  AC,  O1O2  and  BD  intersect 
at  one  point,  denoted  by  O  in  Fig.  11.  Let  a,  (3,  and 
8  denote  the  angles  ZO\OA,  ZOAB  and  ZOO\B, 
respectively. 

It  can  be  shown  that  all  the  circles  {£)}  are  tangential 
to  the  same  boundary  line.  Without  loss  of  generality, 
assume  that  the  tangential  line  is  AC,  i.e.,  di  <  d'i  and 
ri  =  di.  Based  on  simple  geometry,  the  lengths  of  the 
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Fig.  11:  Consider  first  the  downward  part  from  the  grid  to  a  left 
interface  vertex  u~ .  Let  0\  denote  the  location  of  the  last  vertex  on 
the  designated  vertical  line  in  the  grid,  and  O2  the  location  of  v~ . 
The  circles  centered  at  Oi  and  O2  represent  their  maximum  ranges. 
Let  A ,  B  and  C,  D  denote  the  intersecting  points  of  these  two  circles 
with  the  horizontal  lines  at  their  centers.  Let  E±  denote  the  intersection 
between  circle  0 1  and  line  O1O2.  Let  d\  and  denote  the  distance 
between  E\  and  the  two  lines  AC  and  BD,  respectively.  Next  we 
draw  a  circle  with  radius  r\  =  min{di,d^}  centered  at  E\.  Let  E2 
denote  the  intersection  between  circle  E\  and  line  O1O2,  and  a  similar 
circle  centered  at  E2  is  drawn.  This  procedure  is  repeated  to  generate 
a  sequence  of  circles  until  the  last  generated  circle  covers  O2.  This 
sequence  of  circles  {E{\  gives  the  locations  and  the  maximum  ranges 
of  the  vertices  connecting  the  grid  and  v~ .  The  upward  part  from  i/T 
to  the  grid  is  done  with  the  same  procedure  except  starting  from  v~T . 


Fig.  12:  Connecting  the  interface  vertices  with  the  grid. 


line  segments  of  {OEi}  forms  an  equal  ratio  sequence: 

OEi+ 1  =  OEi  -  r *  =  OEi(  1  -  sin  a), 
with  OEi  =  0()\  —  f?i.  We  thus  have 

OEi+ 1  =  (OOi  -  i?i)(l  -  since)*. 

Based  on  the  stopping  condition  of  the  procedure,  the 
number  k  of  circles  is  given  by  the  minimum  index  i 
such  that  OEi+\  <  OO 2.  We  thus  have 
k  =  min{i  £  N  :  OEi+ 1  <  OO2 } 

=  min{i  £  N  :  ( OOi  —  f?i)(l  —  since)*  <  OO2} 
log  — ^iL_ 


=  min{7  eN:i< 


OO1  —  R1 


} 


log(l  —  since)' 

<  log(002/OOi)/ log(l  —  since)  +  1. 


(6) 


Since  AOO2D  and  AOO±B  are  similar  triangles,  the 
ratio  equals  the  ratio  Also  because  —  log(l  — 
x)  >  x  for  0  <  x  <  1,  (6)  can  be  written  as 

k  <  (log Ri  —  log R2)/ since  +  1.  (7) 

Because  0\A0  is  a  triangle  and  /3  >  f ,  the  value  of 


sin  a  can  be  lower  bounded  as  follows: 

.  Ri  .  n  .  Ri  V^2  R\  —  i?2  a/2 


sm  a  = 


■  sin  /3  > 


> 


(8) 


L  <■ 
sin  (3  — 


OOi  00\  2  O1O2  2 

Furthermore,  since  9  =  a  +  (3  >  f3,  the  length  of 

O1O2  =  -jEj  has  an  upper  bound:  O1O2  < 

\[2L.  Hence  (8)  leads  to: 

sin  a  >  (f?i  —  i?2)/(2L).  (9) 

Substituting  (9)  into  (7),  we  have 

k  <  2L(logi?i  -  log  R2) / (7?i  -  R2)  +  1. 
Consider  next  R\  =  i?2-  The  sequence  of  circles  {-E^} 
have  the  same  radius  R\.  Since  O1O2  =  z~Tt  <  2 L,  the 
bound  k  <  +  1  holds.  ■ 

To  satisfy  the  assumption  of  ZCAB  >  j  in  Lemma  4, 
we  set  the  distance  between  the  last  horizontal  line  of 
the  grid  and  the  horizontal  line  of  super  vertices  to  n. 
This  ensures  that  angle  ZB  AO  <  !■  Note  that  in  the 
downward  part  from  the  grid  to  a  left  interface  vertex 


v-  ,  Ri  is  a  constant  and  R2  = 


Hence  the  bound 


on  k  given  in  Lemma  4  can  be  written  as: 


k  < 


< 


2  (log  R\  +  n  log  4  —  log  n) 


-  4-nn/2 


L  +  1 


2(logJ?i  +  nlog4) 


n  +  1, 


i?i  -  1/8 

which  is  in  the  order  of  0(n2).  A  similar  argument  can 
be  made  for  the  upward  part  where  R\  =  4 ~nn  and  R2 
is  a  constant.  The  same  holds  for  R\  =  i?2-  Hence  the 
total  number  of  additional  vertices  to  connect  the  grid  to 
the  interface  vertices  of  a  super  vertex  is  in  the  order  of 
0(n3). 


D.  Reduction  from  MDS  to  TP  in  the  2-D  Disk  Hyper¬ 
graph  H2 

With  H2  constructed,  we  now  establish  the  correctness 
of  the  reduction  from  the  MDS  in  G  to  the  TP  from  v\ 
to  vn+i  in  H2. 

Lemma  5:  Let  ns  =  112  +  1  where  n 2  is  the  total 
number  of  normal  vertices  in  H2-  There  is  a  dominating 
set  with  size  k  in  G  if  and  only  if  there  is  a  path  from  v\ 
to  vn+i  in  H2  with  width  between  kns  and  (k+l)ns  —  1. 

Proof:  The  chosen  value  of  ns  ensures  that  the 
width  of  a  path  from  v±  to  un+i  is  dominated  by  the 
number  of  super  vertices  that  it  covers.  The  correctness 
of  the  reduction  thus  follows  from  the  same  arguments  in 
the  proof  of  Lemma  3  based  on  the  construction  of  H2- 


The  polynomial  nature  of  the  reduction  is  clear  from  the 
construction  of  H2-  We  thus  arrive  at  Theorem  1. 

Appendix  B.  Proof  of  Lemma  1 
Consider  a  TP  problem  from  s  to  t  in  a  k- D  exposed 
disk  hypergraphs  H  =  ( V. ,  E).  We  construct  a  fc-D  UDH 
H'  as  follows.  First,  the  normal  vertex  set  V'  of  //' 
is  given  by  V,  except  that  the  ranges  of  any  v'  £  V' 
equals  max„eyi?t,.  Next,  for  each  vertex  v'  £  V',  we 
place  a  super  vertex  in  (i.e.,  the  exposed  area  of 
the  corresponding  vertex  in  H)  that  contains  \V\  +  1 
normal  vertices  located  sufficiently7  close  to  each  other. 

7The  \V\  +  1  normal  vertices  are  sufficiently  close  such  that  any 
transmission  from  one  of  these  vertices  to  a  vertex  outside  this  super 
vertex  reaches  all  the  \V\  +  1  normal  vertices  in  this  super  vertex. 


II 


The  super  vertices  have  the  same  range  as  the  normal 
vertices  in  V' ,  ensuring  //'  is  a  UDH.  The  reduction  can 
thus  be  seen  by  noticing  that  while  the  enlarged  ranges 
introduce  additional  hyperedges  in  H',  these  hyperedges 
cannot  be  on  a  thinnest  path  due  to  the  fact  that  they  all 
contain  at  least  one  super  vertex. 

Appendix  C.  Proof  of  Lemma  2 

In  this  proof,  we  modify  the  2-D  disk  hypergraph  II-2 
in  the  proof  of  Theorem  1  to  a  2-D  exposed  disk  hy¬ 
pergraph  H3  while  preserving  the  polynomial  reduction. 
Based  on  the  definition,  a  sufficient  condition  for  a  2- 
D  disk  hypergraph  to  be  exposed  is  that  none  of  the 
maximum  range  disks  are  completely  inside  any  other. 
The  vertices  in  H2  for  realizing  the  line  segments  of 
the  grid  satisfy  this  condition.  We  only  need  to  modify 
the  implementations  of  the  crosses  and  around  the  super 
vertices. 

A.  Implementation  of  Crosses 

BC  =  R =  R(J  =  -R  tan  0  ~  0.5543R 

N  BD  =  (1  -  ^rW)R  C8  0.4283K 


Fig.  13:  To  implement  a  directed  cross  shown  in  (a),  we  first 
implement  a  vertex  for  the  blue  line  with  maximum  range  R  at  location 
A  (the  blue  circle)  shown  in  (b).  Next  we  draw  a  perpendicular  bisector 
between  A  and  E  (the  right  intersecting  point  of  the  circle  with 
the  line).  On  this  vertical  line,  we  find  two  points  B  and  C  such 
that  ZBAE  =  Z CAE  =  29°.  At  each  point,  we  put  a  vertex 
for  the  red  line  with  radius  equal  to  the  length  of  BC  (illustrated 
by  the  two  red  circles  in  (b)).  Simple  geometry  calculation  leads  to 
BD  <  BC  <  AB  =  BE.  This  ensures  that  vertices  B  and  C  are 
exposed  yet  cannot  overhear  vertices  located  at  A  and  E.  We  complete 
the  implementation  by  adding  vertices  on  the  vertical  line  BC  and  the 
horizontal  line  AE  (see  (c)  and  (d)).  Note  that  to  preserve  the  exposure 
of  vertices  B  and  C,  the  maximum  ranges  of  vertices  from  point  E  to 
the  right  side  need  to  be  enlarged  gradually  to  the  constant  maximum 
range  of  normal  vertices  on  the  grid  (this  only  requires  a  constant 
number  of  additional  vertices). 

In  the  implementation  of  directed  crosses  in  H2  (see 
Fig.  8),  some  vertices  on  the  line  with  a  lower  level 
index  may  have  an  empty  exposed  area  (see  the  red 
disks  in  Fig.  8  that  are  completely  covered  by  blue 
ones).  To  implement  a  direct  cross  in  a  2-D  exposed 
disk  hypergraph,  the  maximum  ranges  of  vertices  on  the 
line  with  a  lower  level  index  need  to  be  small  enough  to 
preserve  the  direction  of  the  cross  but  also  large  enough 
to  make  the  vertices  exposed.  We  propose  the  scheme 
described  in  Fig.  13. 


B.  Implementation  around  Super  Vertices 

In  the  previous  implementation  around  a  super  vertex 
Vj,  all  the  vertices  are  exposed  except  the  anchor  vertices 
{u~,u~l}  and  the  bridging  vertices  {pik}-  However, 
we  notice  that  these  vertices  would  all  be  exposed  if 
there  were  no  interface  vertices.  Our  solution  is  thus  to 
move  all  the  interface  vertices  away  from  their  original 
positions  by  a  constant  distance  and  add  a  constant 
number  of  vertices  to  connect  each  new  interface  vertex 
to  the  bridging  vertex  or  the  anchor  vertex  on  the  right 
side.  A  detailed  implementation  is  shown  in  Fig.  14. 


Fig.  14:  An  interface  vertex  on  the  left  side  is  replaced  by  three 
vertices  with  maximum  ranges  r,  fr  and  f  respectively.  These  three 
vertices  are  located  on  a  vertical  line  to  the  left  side  of  the  original 
location  of  the  interface  vertex  with  a  distance  of  r.  An  interface  vertex 
on  the  right  side  is  replaced  by  two  vertices  with  maximum  range 
located  on  a  vertical  line  to  the  right  side  of  the  original  location  of 
the  interface  vertex  with  distance  J.  Under  this  implementation,  the 
exposed  areas  of  the  anchor  and  bridging  vertices  are  right  above  the 
point  where  they  are  tangential  with  the  horizontal  line  of  the  super 
vertices  (as  illustrated  by  the  arrows). 

Appendix  D.  Proof  of  Theorem  4 

Consider  an  MDS  problem  in  a  graph  G  with  a 
maximum  degree  of  3.  We  first  follow  the  first  two  steps 
in  the  proof  of  Theorem  1  to  build  the  grid  representation 
of  hypergraph  Hi.  Note  that  due  to  the  unit  range  of  all 
vertices,  we  set  the  size  of  the  grid  to  a  constant  greater 
than  1  (say,  5)  to  avoid  unwanted  overhearing.  Next, 
we  implement  this  representation  in  a  3-D  UDG  while 
preserving  the  reduction.  Any  line  segment  of  hyper¬ 
edges  in  Hi  is  replaced  by  a  sequence  of  unit  disks,  one 
just  touching  the  another.  Any  cross  between  two  line 
segments  can  be  easily  implemented  by  using  the  third 
dimension,  as  shown  in  Fig.  15.  In  this  implementation, 
there  is  no  overhearing  between  vertices  on  these  two 
line  segments  at  all.  Since  G  has  a  maximum  degree 
of  3,  there  are  at  most  4  hyperedges  passing  through 
a  super  vertex.  It  can  be  easily  implemented  without 
any  unwanted  overhearing  (see  Fig.  16).  To  prevent  the 
super  vertices  from  relaying  messages,  we  place  a  mega 
vertex  besides  each  super  vertex.  This  mega  vertex  is 
only  within  the  range  of  this  super  vertex  and  contains 
more  normal  vertices  than  the  total  number  of  normal 
vertices  in  the  reduced  graph  (including  the  normal 
vertices  contained  in  all  the  super  vertices  but  not  those 
in  other  mega  vertices).  In  this  way,  a  path  via  any  super 
vertex  covers  at  least  one  mega  vertex,  thus  cannot  be 
the  thinnest  path.  Fig.  16  illustrates  the  implementation 
around  a  super  vertex8.  The  correctness  of  the  reduction 

8We  can  consider  reduction  from  MDS  in  graphs  with  a  maximum 
degree  up  to  9.  In  this  case,  there  are  at  most  10  incoming  hyperedges. 
Along  with  the  mega  vertex,  they  can  be  packed  around  a  super  vertex 
without  overhearing. 


follows  from  the  same  arguments  as  in  the  proof  of 
Lemma  3. 


Fig.  15:  An  implementation  of  a  cross  in  3-D  UDG. 


Appendix  E.  Proof  of  Theorem  5 

We  first  show  that  as  long  as  there  exists  a  path  from 
s  to  t,  there  exists  a  path  from  s  to  u/  that  traverses  only 
the  sub-hypergraph  H' .  This  can  be  shown  by  noticing 
that  ui  must  hear  the  message  from  s  before  U[-\  and 
any  vertex  to  the  right  of  ui- 1.  This  is  due  to  the 
monotonicity  of  wireless  broadcast  and  the  definition  of 
predecessor.  Consequently,  there  must  exist  a  path  from 
s  to  ui  in  //'.  Since  V  is  covered  by  the  hyperedge 
leading  from  ui  to  iq_i  in  L\,  the  concatenation  of  L\ 
with  any  path  to  m  in  H'  covers  the  same  set  of  vertices. 
Specifically,  the  cover  of  the  path  returned  by  NBI  is  the 
set  of  vertices  located  between  (and  including)  ui  and  t. 
Since  any  path  from  s  to  t  covers  this  set  of  vertices, 
the  correctness  of  the  algorithm  is  established. 

Next,  we  prove  the  property  of  Al*  under  the  disk 
propagation  model.  We  first  state  the  following  lemma 
that  follows  directly  from  triangle  inequality. 

Lemma  6:  Let  D\  and  I) 2  denote  two  closed  balls 
in  IR/y  with  radii  r\  and  r^,  respectively.  Let  a  denote 
the  distance  between  the  centers  of  D\  and  D^.  If  0  < 
a  <  |ri  —  r2 1 ,  then  D2  C  D\. 

Based  on  Lemma  6,  for  any  vertex  v  between  ui  and 
ui- 1,  we  have  Dv>Rv  C  AUl^Ul Therefore  AL,  = 

ALi  =  ULi  Duk,d(uk,uk _i)  (let  u0  =  t ).  Next,  consider 
an  arbitrary  path  L  from  s  to  t.  We  show  that  for  any  Uk 
(k  =  1, . . . ,/),  £>Ufc,d(ufc, Uk_!)  C  Al.  Specifically,  since 
Uk- 1  must  first  hear  the  message  from  Uk  or  a  vertex  to 
the  left  of  Uk,  DUk,d{uk,uk_1)  is  a  subset  of  the  covered 
area  of  this  hop  in  L  based  on  Lemma  6.  This  completes 
the  proof. 

Appendix  F.  Proof  of  Theorem  7 
A.  For  General  Directed  Flypergraphs 

Let  L\  denote  the  path  from  s  to  f  provided  by 
SPBA  and  Lopt  =  {ei,  e-z,  ■  ■  ■ ,  ek}  the  thinnest  path. 
If  multiple  thinnest  paths  exist,  let  Lopt  be  the  one  with 
the  minimum  number  of  hyperedges.  Let  C(L)  denote 
the  length  (i.e.,  the  sum  of  hyperedge  weights)  of  L. 

Since  each  vertex  covered  in  L  \  (except  the  source  s) 
contributes  to  the  weight  of  at  least  one  hyperedge  in 
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L 1,  the  width  W{L\ )  is  no  larger  than  the  length  of  this 
path  plus  one.  Also  because  L\  is  the  shortest  path,  its 
length  is  no  larger  than  the  length  of  Lopt.  We  thus  have 
W(Li)  <  C(Li)  +  1  <  C(Lopt)  +  1.  (10) 

We  then  obtain  the  approximation  ratio  by  deriving  an 
upper  bound  of  C(Lopt )  as  a  function  of  W(Lopt). 

Note  that  the  destination  set  Te  of  hyperedge  e,  on 
Lopt  cannot  contain  k  —  i  +  1  vertices:  its  own  source 
vertex  sei  and  vertices  in  {sei+2 ,  sei+3 , . . . ,  sek ,  t }.  The 
later  holds  because  otherwise  Lopt  is  not  the  the  thinnest 

path  with  minimum  number  of  hyperedges.  We  thus  have 

k 

C(Lopt)  <  J2(W(Lopt)  -{k-i  +  1)) 

=  kW{Lopt)  -  k(k  +  l)/2 
<  W(Lopt)(W(Lopt)  -  l)/2,  (11) 

where  (11)  comes  from  k  <  W(Lopt )  —  1.  Substituting 
(11)  into  (10),  we  have 

W{Li)  <  W{Lopt)(W(Lopt)  -  l)/2  +  1 

<  W2{Lopt)/ 2,  (12) 

where  (12)  holds  since  W(Lopt)  >  2. 

Based  on  (12),  if  W(Lopt)  <  \/2n,  then  W(L\)  < 
\ W2  (Lopt)  <  Z^W{Lopt).  Otherwise,  we  have 
W(L\)  <  n  <  ,y^W(Lopt).  In  summary,  SPBA 

provides  a  approximation. 

B.  For  Ring  Flypergraphs 

Since  a  ring  hypergraph  is  a  special  directed  hyper¬ 
graph,  all  the  analysis  in  the  previous  subsection  applies. 
Specifically,  inequality  (10)  holds.  The  problem  then 
remains  in  obtaining  a  tighter  upper  bound  of  C(Lopt) 
based  on  the  geometrical  properties  of  ring  hypergraphs. 

First,  note  that  the  length  of  a  hyperpath  L  equals  the 
sum  of  the  number  of  times  each  vertex  is  reached.  Let 
Ev  denote  the  set  of  hyperedges  on  Lopt  that  include  v 
in  their  destination  sets,  i.e., 

Ev={e  &  Lopt  '■  v  £  Te}- 

Now  we  construct  a  subset  E'v  of  Ev  by  iteratively  re¬ 
moving  one  from  any  pair  of  hyperedges  whose  positions 
in  Lopt  are  adjacent  until  no  such  pair  exists.  Because  at 
most  half  of  the  hyperedges  are  removed  from  Ev,  the 
size  of  E'v  is  at  least  half  of  the  size  of  Ev,  in  another 
word  |^|  <  2\E'V\. 

Let  Rmax  and  Rmin  denote  the  largest  maximum 
range  and  the  smallest  minimum  range  among  all  ver¬ 
tices  in  the  given  ring  hypergraph  Hr,  respectively.  Let 
R'rnin  be  the  larger  one  between  Rmin  and  the  smallest 
distance  between  any  two  vertices  in  Flr.  Based  on 
the  construction  of  E'v,  the  set  of  source  vertices  of 
hyperedges  in  E'v  satisfies  two  properties.  First,  based  on 
the  definition  of  ring  hypergraphs,  the  distance  between 
any  source  vertex  in  the  set  and  v  is  no  larger  than  the 
maximum  range  of  this  vertex  and  hence  no  larger  than 
Rmax ■  Second,  the  distances  between  any  two  source 
vertices  in  the  set  are  larger  than  Rmin  and  hence 
R'min-  Otherwise  the  two  hyperedges  rooted  at  these  two 
vertices  can  reach  the  source  vertex  of  each  other  and 
hence  they  are  adjacent  in  Lopt  (recall  that  e,  £  Lopt 
cannot  reach  any  vertex  in  {sei+2,  sei+3, . . . ,  sefe}). 


Mega  vertex 

Fig.  16:  Implementation  around  the  super  vertices  in  UDG. 
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Given  these  two  properties,  the  size  of  E'v  thus  is 
upper  bounded  by  the  maximum  number  of  points  in  the 
Euclidean  space  that  are  at  most  Rmax  away  from  v  and 
at  least  R'min  apart  from  each  other.  This  is  equivalent 
to  a  sphere  packing  problem  of  arranging  the  maximum 
number  of  small  spheres  with  radius  R'min/2  inside  a 
large  sphere  with  radius  Rmax  +  R'min/2.  An  upper 
bound  of  this  packing  problem  is  the  ratio  between  the 
volumes  of  the  large  and  small  spheres.  We  thus  have: 


\K\< 


( Rmax  +  R'min/^)0 


=  (1  +  2a)d, 


(R'miJVd 

where  a  =  Rmax / R'min.  Recall  that  \EV\  <  2\E'V\.  Note 
that  the  destination  t  can  only  be  reached  by  the  last 
hyperedge  e*,  and  hence  \Et\  =  I {efe} I  =  1-  We  thus 
have 

C(Lopt)=  J2  (13) 

V& Lopt\{t} 


<  2(1  +  2 a)d{W(Lopt)  -  1)  +  1 

<  2(1  +  2a)dW{Lopt)  -  1.  (14) 

Substituting  (14)  into  (10).  we  have 

W(L i)  <  C(Lopt)  +  1  <  2(1  +  2 a)dW{Lopt),  (15) 
i.e.,  SPBA  provides  a  2(1  +  2a)d-approximation  for  TP 
in  d-D  ring  hypergrpahs. 


C.  Asymptotic  tightness 

We  now  prove  that  y^-ratio  is  asymptotically  tight 
even  for  2-D  disk  hypergrpahs.  The  proof  has  two  steps. 
First,  we  construct  a  directed  hypergraph  H  for  which 
the  worst  case  ratio  is  asymptotically  reached.  Next,  we 
show  a  2-D  disk  implementation  of  H. 

Consider  the  the  following  hypergraph  H  illustrated  in 
Fig.  17  with  k  red  vertices  v\ , . . .  ,Vk  and  k'  blue  vertices 
Mi , . . . ,  Mfc'  along  with  the  source  s  and  the  destination 
t.  Each  red  vertex  vt  has  one  outgoing  hyperedge  e  with 
Te  =  {mi,  . . . ,  Vi-i,  Mj+i}  (let  Vk+  i  denote  t).  Each  blue 
Mj  has  one  outgoing  hyperedge  e  with  Te  =  {Mi+i}  (let 
Mfc'+i  denote  t).  Finally,  we  add  two  hyperedges  that 
connect  source  s  to  v\  and  mi  respectively. 


Fig.  17:  A  worst  case  scenario  for  SPBA. 

Let  k!  =  k(k  +  l)/2  +  1.  Since  the  shortest  path 
traverses  through  the  blue  hyperedges  while  the  thinnest 
path  through  the  red  ones,  the  approximation  ratio  is 
given  by: 

7(fc)  =  (k2  +  k  +  2)/{2k  +  4).  (16) 

Note  that  the  total  number  of  vertices  is 

n  =  k  +  k'  +  2  =  k  +  k(k  +  l)/2  +  1. 

When  n  is  large,  k  ~  \[2 n  and  a(k)  ~  -|  ~  \J~%- 
Next,  we  implement  the  above  hypergraph  under  a  2-D 
disk  model  as  illustrated  in  Fig.  18.  The  red  vertices  are 
located  on  a  straight  line  with  RVl  =  RV2  =  1,  RVi  = 
2l~ 2  for  i  >  2.  The  source  vertex  s  is  located  on  the 
line  to  the  left  of  mi ,  and  both  its  maximum  range  and  its 
distance  to  v\  equals  to  RVk .  The  terminal  vertex  t  has  a 


maximum  range  of  0  and  is  located  to  the  right  of  Vk  with 
a  distance  of  RVk .  The  maximum  range  of  a  blue  vertex 
Ui  is  RVk  —  ie  where  e  is  a  small  positive  value  to  prevent 
Ui-i  from  overhearing  messages  transmitted  by  Ui.  And 
the  blue  vertices  are  located  on  a  route  from  s  to  t  that 
contains  two  vertical  line  segments  of  length  (1  +  l)RVk 
and  a  horizontal  one  of  length  3 RVk ,  as  demonstrated  by 
the  blue  dashed  lines  in  Fig.  18.  The  positive  parameter 
l  is  used  to  prevent  a  blue  vertex  from  overhearing  the 
last  red  vertex  Vk ■  In  the  asymptotic  regime  with  large  k, 
l  can  be  set  sufficiently  large  so  that  the  k'  blue  vertices 
can  be  implemented  along  the  depicted  route  from  s  to 
t. 


SPBA. 


Appendix  G.  Proof  of  Theorem  8 

Let  L2  denote  the  path  in  hypergraph  H  from  s  to 
t  given  by  the  TSBA  algortihm  and  Lopt  the  thinnest 
path.  Let  Li(v)  and  L2(v)  denote  the  paths  from  s  to 
a  vertex  v  given  by  SPBA  and  TSBA,  respectively.  The 
following  lemma  establishes  a  property  of  L2(y). 

Lemma  7:  For  any  hyperedge  e  in  H,  we  have,  Vm  € 

Te, 

W(L2(v ))  <  \L2(se)  U  Te\. 

Proof:  Lemma  7  follows  directly  from  the  tree 
structure  of  TSBA.  ■ 


A.  For  General  Directed  Hypergraphs 

Let  Lopt  =  {ei, . . . ,  e*,}  denote  the  thinnest  path.  For 
ease  of  presentation,  let  the  sequence  of  source  vertices 
sei , . . . ,  sek  and  the  final  destination  t  be  denoted  as 
Mi, ... ,  Vk+ 1-  Let  U  =  Based  on  Lemma  7,  we 

have,  for  all  i  =  1, . . . ,  k, 

W(L2(vi+1))  <  \L2(yi)  U  Tei  | 

<  W(L2(vi))  +  \Tei\{Vl, . . . ,  m1+i}|  +  1 
=  W(L2{Vi))  +  \Tei\U\  +  l,  (17) 

where  (17)  holds  since  Tei  does  not  contain  vertices  in 
{sei+2,  •  ■ . ,  sek,t}.  Summing  (17)  over  i,  and  noticing 

that  ^(/^(fi))  =  1  and  L2(vk+ 1)  =  L2,  we  have: 

k 

W(L2)<k  +  l  +  J2\TeAU\- 

Next,  since 

W(Lopt)  =  k+  1  +  luJU  (Tei\U)\, 
we  can  upper  bound  \Tei\U\  by  W(Lopt)  —  k  —  1  for 
any  i  =  1 , ,k.  Thus 

W{L2)  <  k  +  1  +  k(W(Lopt)  -k-  1). 
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The  right  side  of  this  inequality  is  a  quadratic  function 
of  k  with  the  maximum  at  k  =  W(Lopt)/ 2.  We  thus 
have 

W(L2)  <l  +  W2{Lopt)/4. 

If  W(Lopt)  <  2V^1,  the  approximation  ratio  7  is 
given  by 

7  =  W(Lopt)/ 4  +  1  /W(Lopt)  <  n/ 2jW=l.  (18) 
The  inequality  holds  because  the  function  |  +  is  an 
increasing  function  for  x  >  2. 

If  W(Lopt)  >  2 y/n  —  1,  we  have 

7  <  n/W(Lopt )  <  n/{ 2V/^I).  (19) 

This  completes  the  proof  for  case  of  general  directed 
hypergraphs. 


B.  For  Ring  Hypergraplis 

Let  Li(t)  =  {ei,  e2, . . . ,  &k]  be  the  shortest  path  from 
s  to  t.  Let  Vi  denote  the  source  vertex  of  e,  and  Vk+i  =  t. 
We  prove,  through  induction,  the  following  inequality  for 
all  i  =  1, . . . ,  k  +  1: 

W(L2(Vi))  <£(Li(Vi))  +  l.  (20) 

When  i  =  1,  (20)  holds  since 

W(L2(v1))  =  1,C(L1(v1))=0. 

Now  assume  that  (20)  holds  for  i  —  1,  i.e., 
W (L2(vi-i))  <  C(Li(vi-i ))  +  1.  Based  on  Lemma  7 
and  this  induction  assumption,  we  have: 

W(L2(Vi))  <  lUivi-JUTe^l 

<  W (L2(Vi_i))  +  \Tei_1\ 

<  +  1  +  |Tei_i| 

=  C(L\(vi))  +  1. 

This  completes  the  induction.  Considering  Vk+i  =  t,  we 
have 

W(L2)<C(L1(t))  + 1.  (21) 

From  (15),  (10),  and  (21),  we  have  W(L2(t))  < 
2(1  +  2 a)dW{Lopt),  i.e.,  TSBA  provides  a  2(1  +  2a)2 
approximation  for  ring  hypergrpahs. 

C.  Asymptotic  tightness 

We  first  construct  a  directed  hypergraph  H  as  illus¬ 
trated  in  Fig.  19.  The  vertex  set  of  H  consists  of  two 
types  of  vertices:  k  +  1  normal  vertices  vo,  .  ■  ■ ,  Vk-i  and 
t,  and  k  super  vertices  ui, . . . ,  Uk-i,u,  each  containing 
k  —  1  normal  vertices.  Rooted  at  each  normal  vertex  17 
are  two  hyperedges  e,;+i  and  e'i+1.  Hyperedge  e,  has 
destination  vertices  Tei  =  {17,17}  and  hyperedge  e'  has 
destination  vertices  Te>  =  {17,11}. 

It  is  easy  to  see  that  the  thinnest  path  from  vo  to  t 
is  Lopt  =  {e},  e'2,  ■  ■  ■ ,  e'fe,  e'k}  with  width  2k.  However, 
TSBA  returns  the  path  Ld  =  {ei,  e2: . . . ,  e&}  with  width 
k2  +  1  in  the  worst  case9.  The  approximation  ratio  is 
W(L)  k2  + 1  n 

7  =  W(Lopt)  =  2k  =  V^T' 

Given  the  similarity  between  H  and  the  hypergraph 
Hi  constructed  in  the  proof  of  Theorem  1,  we  can 
follow  the  same  approach  given  in  Appendix  A  to 

9Note  that  Vi  can  update  1  through  both  and  with 

the  same  width.  Since  the  order  of  hyperedges  used  in  the  update  is 
arbitrary,  \  could  be  used  to  update  174.1  for  all  I  <  i  <  k  ! 
in  the  worst  case. 


implement  H  under  a  2-D  disk  model.  However,  this 
implementation  requires  additional  vertices  (referred  to 
as  auxiliary  vertices)  that  may  render  our  previous 
approximation  analysis  invalid.  To  maintain  the  ratio, 
each  original  vertex  (including  the  vertices  in  a  super 
vertex)  in  H  is  replaced  with  c  vertices  (clustered 
together)  in  its  2-D  disk  implementation,  where  c  is  the 
number  of  auxiliary  vertices  introduced  by  the  imple¬ 
mentation.  In  this  case,  TSBA  returns  a  path  that  covers 
{111,  ■  ■  ■ ,  itfe_  1,17), . . . ,  17c}  along  with  a  set  of  auxiliary 
vertices.  The  thinnest  path  covers  {it, no ,...,Vk}  and 
another  set  of  auxiliary  vertices.  The  approximation  ratio 
in  this  2-D  disk  hypergraph  is  given  by 

_  (fc2  +  l)c  +  c/  _  k2  +  1  +  ^ 

7  2  kc  +  c"  2k  +  — 

C 

where  d  and  c"  denote  the  number  of  auxiliary  vertices 

covered  by  the  path  returned  by  TSBA  and  the  thinnest 

/  // 

path.  Since  —  <  1  and  —  <  1,  when  n  is  large,  we 
have  7  ~  2y/n-v  7e'’  t^le  aPProximation  ratio  is 

asymptotically  tight  in  2-D  disk  hypergraphs. 
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